† Corresponding author. E-mail:
Project supported by the National Science and Technology Pillar Program, China (Grant No. 2015BAK17B06), the Earthquake Industry Special Science Research Foundation Project, China (Grant No. 201508026-02), the Natural Science Foundation of Heilongjiang Province, China (Grant No. A201310), and the Scientific Research Starting Foundation for Post Doctorate of Heilongjiang Province, China (Grant No. LBHQ13040).
Twin gold crystal nanowires, whose loading direction is parallel to the twin boundary orientation, are simulated. We calculate the nanowires under tensile or compressive loads, different length nanowires, and different twin boundary nanowires respectively. The Young modulus of nanowires under compressive load is about twice that under tensile load. The compressive properties of twin gold nanowires are superior to their tensile properties. For different length nanowires, there is a critical value of length with respect to the mechanical properties. When the length of nanowire is greater than the critical value, its mechanical properties are sensitive to length. The twin boundary spacing hardly affects the mechanical properties.
In recent years, metal nanomaterials, whose properties are different from macroscopic metal materials, have received a great deal of attention. The experiment of metal nanowires is uncontrollable, but by the molecular dynamics (MD) method their deformation can be observed under a certain condition.
Gleiter[1,2] investigated the deformation of nano-structure. Janusz et al.[3] used the molecular dynamics simulation method to calculate thermodynamic properties of the transition metal with face-centered cubic structure. Park and Zimmerman[4] discussed the mechanical properties of the gold nanowires in the yield condition, and summarized the plastic characteristics of gold nanowires by molecular dynamics simulation.
The current researches focus on the twin structure, whose load direction is vertical to the twin boundary orientation. The point defects in Au twin boundaries,[5] the fracture behavior and the ductility of the twinned Cu nanowires,[6] and the ductility of the magnesium nanowire[7] were studied under loads in the above direction.
In the present paper, we discuss the properties of gold twin crystal nanowires, whose directions are all parallel to the twin boundary orientation, under different conditions.
The rest of the paper is organized as follows. In Section 2, we discuss the simulation model. The results will be discussed in Section 3. Finally, some conclusions are drawn from the present study in Section 4.
Some twin gold crystal models are given by the geometric construction method (Fig.
The molecular dynamics simulation software NanoMD,[8] which is developed by Zhao’s group in Nanjing University, are modified[9,10] for simulation. In this paper we adopt Johnson’s EAM potential function[11] to describe the interaction between each two gold atoms. Cell list and Verlet list are combined with integration step δt = 2.94 fs = 2.94× 10−15 s. At the beginning, atom initial velocity, which fits the Maxwell–Boltzmann distribution,[12] is generated randomly. The temperature of the model system keeps up 300 K by the Nosé-Hoover[13–16] method. Through 2×104 relaxation steps, the system achieves a thermal equilibrium state. The loads is exerted on the [110] orientation in the z-axial direction, whose stress is obtained by[17]
The nanowires under tensile or compressive loads and different lengths are simulated respectively. Their mechanical properties are different under different conditions.
The models deform at 0.129 m/s tensile or compressive velocity rate in the z-axial direction on the basis of the above mentioned models.
The average energies of models under different load styles are different (Fig.
The average potential energies in different steps under tensile and compressive loads are shown in Fig.
Figure
Table
In the paper, we discuss the micro-deformations by using colored atom balls, where blue atoms are the twin boundary or dislocation atoms and yellow atoms are stable.
Relaxations under tensile and compressive loads are similar due to the fact that models are not applied loads. In the process of compression, compressive instability leads to the fact that dislocations appear, which come into slip planes. The models will be broken in the form of pileup along slip planes. Nevertheless, there are different statuses in the models under tensile loads. Slip planes break the twin crystal structure, around which the non-crystallizing of the atoms is acute and the necking phenomenon turns up until it is broken. The directions of dislocations are both 45 angle and twin boundaries prevent the dislocation from developing.
On the basis of the above mentioned models, different length twin gold nanowires are simulated. Different mechanical properties are obtained.
The energy increases linearly during the relaxation. The system’s average kinetic energy is dependent only on temperature (Eq. (
Although the lengths of gold nanowires are different, the variation trends of curves are similar (Fig.
Figure
Also, in the progress of tension it is verified that the dislocations destroy the prevention of a twin boundary and keep developing along the length direction.[18]
The twin boundary spacing has an effect on the mechanical properties when the tensile direction is vertical to the twin boundary orientation.[19] There appear different properties when the loading direction is parallel. The other two twin boundary models are simulated (Fig.
Figure
In this work, we investigate the mechanical properties of twin gold crystal nanowires by the molecular dynamics method through controlling variables. The nanowires under tensile or compressive loads, different length nanowires, and different twin boundary nanowires are simulated respectively. The micro-deformation is the progress of dislocation developing by breaking twin boundaries.
When the models are under different style loads, the mechanical properties are different. The Yong modulus under a compressive load is about twice that under a tensile load. The compressive properties of twin gold nanowires are superior to its tensile properties.
The length of nanowires also affects the mechanical properties. When the length of nanowires is short: L = 10a, 13a, 16a, the nanowires length changes hardly affect the mechanical properties of twin gold nanowires; when the length turns long with the increase of the nanowire length, system yield stress and Young modulus decrease and the yield strain increases. The longer the nanowires, the less the influence of non-crystallization is.
When the loading direction is parallel to the twin boundary, the twin boundary spacing hardly affects the mechanical properties.
1 | |
2 | |
3 | |
4 | |
5 | |
6 | |
7 | |
8 | |
9 | |
10 | |
11 | |
12 | |
13 | |
14 | |
15 | |
16 | |
17 | |
18 | |
19 |